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ABSTRACT 


The  electromagnetic  field  generated  by  a progressive  ocean 
wave  in  a horizontally  stratified  ocean  is  a sum  of  a transverse  elec- 
tric type  field,  a transverse  magnetic  type  field,  and  an  electrostatic 
type  field.  Sea  water  velocity  components  in  a vertical  plane  contain- 
ing the  direction  of  wave  propagation  generate  the  transverse  electric 
part  of  the  field,  and  the  velocity  component  normal  to  the  plane  gen- 
erates the  transverse  magnetic  part  of  the  field,  which  vanishes  above 
an  ocean  surface.  The  electrostatic  part  of  the  field  results  from  sur- 
face charges  that  halt  vertical  electric  currents  at  an  ocean  surface. 

Gradients  of  magnetiv  fields  induced  at  the  surface  by  both 
surface  and  internal  waves  provide  sensible  signals  for  recently 
developed  magnetic  gradiometers  based  on  the  Josephson  effect 
provided  surface  speeds  exceed  1 cm/sec  or  so.  Vertically  spaced 
measurements  of  noise  spectra  of  magnetic  field  gradients  above  an 
ocean  surface  offer  a unique  and  promising  means  of  obtaining  a 
measure  of  surface  and  internal  wave  spectra  for  three  reasons: 

(1)  magnetic  field  strengths  above  the  surface  are  proportional 
to  a weighted  average  of  sea  water  speed  over  an  ocean  depth,  which 
provides  a response  depending  on  the  mode  structure  of  internal 
waves,  (2)  field  strengths  decrease  exponentially  with  height  above 
the  surface  as  e , where  k denotes  wave  number,  which  provides 
a means  of  wave  number  discrimination  and,  together  with  frequency 
discrimination,  offers  means  of  resolving  internal  wave  spectra  mode 
by  mode,  and  (3)  noise  spectra  of  magnetic  field  gradients  are  effec- 
tively cross-spectra  and  so  provide  directional  information  on  wave 
spectra  from  an  effective  point  measurement. 
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INTRODUCTION 


As  Faraday  [1832]  first  recognized,  sea  water  moving 
across  the  Earth's  steady  magnetic  field  induces  eddy  currents  in 
a sea  that  generate  an  electromagnetic  field  characteristic  of  the 
velocity  field  in  the  sea.  First  investigations  of  the  effect  were 
concerned  with  measurement  of  electric  fields  induced  by  steady 
ocean  currents.  Efforts  culminated  in  development  of  a so-called 
geomagnetic  electrokinetographic  or  GEK  method  [von  Arx,  1950] 
of  inferring  velocities  of  broad,  steady,  ocean  currents  from 
measurements  of  horizontal  differences  in  electric  potential  at  an 
ocean  surface. 

Longuet-Higgins,  et  al.  [1954]  develop  a sound  basis 
for  interpreting  GEK  measurements  and  consider  electric  fields 
and  currents  produced  by  surface  waves  as  well  as  steady  currents. 
Investigation  of  magnetic  fields  generated  by  either  steady  ocean 
currents  or  ocean  waves,  however,  has  been  largely  neglected  until 
recently.  Larsen  [1968]  and  Sanford  [1971]  among  others  investigate 
electric  and  magnetic  fields  generated  by  deep  sea  tides  and  other 
large scaJe,  quasistatic,  ocean  currents.  Crews  and  Futterman  [1962] 
are  the  first  to  investigate  magnetic  fields  produced  by  surface  waves. 
They  neglect  magnetic  induction  effects  and,  by  using  the  integral 
form  of  Ampere’s  law,  determine  magnetic  fields  produced  above  the 
surface  by  waves  on  the  surface  of  a deep  ocean.  Warburton  and 
Caminiti  [1964]  extend  the  analysis  to  determine  magnetic  fields 
produced  below  the  surface  of  a deep  ocean,  and  Groskaya,  et  al. 
[1972]  determine  magnetic  fields  produced  by  surface  waves  in 
shallow  seas.  Weaver  [1965]  is  the  first  to  account  for  magnetic 


induction  effects  in  determining  magnetic  fields  induced  by  surface 
waves  on  a deep  ocean  by  using  differential  forms  of  both  Faraday's 
and  Ampere's  laws.  Beal  and  Weaver  [1970]  also  investigate  mag- 
netic fields  produced  by  internal  waves  propagating  along  a sharp 
thermocline.  Nonetheless,  their  investigations  are  restricted  to 
irrotational  velocity  fields,  which  is  a suitable  approximation  for 
surface  waves  but  unnecessarily  restricts  results  for  internal  waves 
and  precludes  application  to  waves  influenced  by  inertial  or  Coriolis 
forces. 

A primary  purpose  of  this  paper  is  to  generalize  and 
extend  previous  work  in  order  to  remove  restrictive  assumptions 
imposed  on  sea  water  velocity  fields.  We  presume  only  that  the  diver- 
gence of  velocity  vanishes  and  so  only  exclude  flow  fields  produced  by 
sound  waves.  By  using  Faraday's  and  Ampere's  laws  to  derive  a vr  .tor 
potential  function  for  the  induced  electromagnetic  field,  we  demon- 
strate that  the  electromagnetic  field  generated  by  a progressive  ocean 
wave  in  a horizontally  stratified  ocean  is  a sum  of  a transverse 
electric  type  field,  a transverse  magnetic  type  field,  and  an  electro- 
static field.  Velocity  components  in  a vertical  plane  containing  the 
direction  of  wave  propagation  generate  the  transverse  electric 
part  of  the  field,  and  the  velocity  component  normal  to  the  plane 
generates  the  transverse  magnetic  part  of  the  field.  The  electro- 
static part  of  the  field  results  from  surface  charges  that  halt  vertical 
electric  currents  at  an  ocean  surface. 

At  wave  frequencies  appreciably  greater  than  tidal  frequen- 
cies, sea  water  velocities  are  elliptically  polarized  in  a vertical 
plane  containing  the  direction  of  wave  propagation.  Consequently, 
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surface  waves  and  high  frequency  internal  waves  generate  transverse 
electric  type  fields  alone,  accompanied  by  an  electrostatic  field.  Both 
transverse  electric  and  transverse  magnetic  type  fields,  together 
with  accompanying  electrostatic  fields,  are  generated  by  ocean  waves 
with  frequencies  of  the  order  of  tidal  frequencies.  Transverse  magnetic 
type  fields,  however,  vanish  above  an  ocean  surface  and  so  are  confined 
below  the  surface,  but  electrostatic  and  transverse  electric  type  fields 
are  present  both  above  and  below  the  surface. 

A second  purpose  of  this  paper  is  to  put  forth  and  examine 
means  of  obtaining  a measure  of  ocean  wave  spectra,  for  both  surface 
and  internal  waves,  by  measuring  noise  spectra  of  magnetic  field 
gradients  above  an  ocean  surface.*  The  magnetic  field  generated  by 
moving  sea  water  is  a linear  transformation  of  the  sea  water  velocity 
field  and  so  is  effectively  the  response  of  a linear  network  — the  Earth's 
magnetic  field.  Consequently,  a correspondence  between  noise  spectra 
of  magnetic  fields  and,  hence,  magnetic  field  gradients  induced  by 
ocean  waves  and  velocity  spectra  of  corresponding  wave  fields  is 
provided  by  a transfer  function  that  is,  in  effect,  the  magnetic  field 
induced  by  a progressive  ocean  wave  of  unit  magnitude. 

At  present,  measurements  of  magnetic  fields  over  oceans 
are  largely  limited  to  mapping  large  scale  irregular  features  of  the 
Earth's  field  [Bullard  and  Mason,  1963,  and  Heirtzler,  1970]  because 
of  limitations  on  sensitivity  and  difficulties  of  separating  spatial  and 
temporal  fluctuations  produced  by  ocean  waves,  ionospheric  sources, 
and  geologic  sources.  Nonetheless,  recent  measurements  of  magnetic 


♦Measuring  spatial  gradients  in  magnetic  fields  rather  than  field 
magnitudes  suppresses  noise  produced  by  ionospheric  and  other 
remote  sources. 
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fields  induced  by  surface  waves  using  so-called  total  field  or  alkali 

vapor  magnetometers  on  fixed  and  floating  platforms  are  reported 

by  Maclure,  et  al.  [1964],  Kozlov,  et  al.  [1971],  Kravtsov,  et  al. 

[1971],  and  Grosk*ya,  et  al.  [1972],  Baker  and  Graefe,  [1968]  report 

limited  success  in  measuring  surface  wave  spectra  using  a rubidium 

vapor  magnetometer  towed  by  aircraft.  The  advent  of  superconducting 

magnetometers  based  on  the  Josephson  effect  [Clarke,  1974  and  Goodman, 

et  al.  1973],  however,  offers  unprecedented  sensitivities  (10~^  gamma/ 

/hertz)  that  effectively  allow  measurement  of  field  gradients  at  a point 

with  sensitivities  of  the  order  of  10"*  (gamma/ m)// hertz,  with  the 

-9  2 

definition  one  gamma  equals  10  webers/m  . 

In  order  to  examine  potential  utility  of  superconducting 
magnetic  gradiometers  for  measuring  ocean  wave  spectra,  we  first 
determine  magnetic  fields  induced  by  progressive  ocean  waves  and 
so  specify  transfer  functions  relating  noise  spectra  of  magnetic  fields 
and  corresponding  ocean  wave  spectra.  We  find  that  magnetic  fields 
induced  above  the  surface  by  progressive  waves  are  circularly 
polarized  in  a vertical  plane  containing  the  direction  of  wave  propa- 
gation and  that  their  magnitudes  decrease  exponentially  with  height  above  the 
surface  as  e , where  k denotes  wave  number  of  a progressive  wave. 
Magnitudes  at  the  surface  are  proportional  to  a weighted  average  of 
sea  water  speed  over  an  ocean  depth.  Magnetic  field  gradients  induced 
at  the  surface  by  both  surface  and  internal  waves  provide  sensible 
signals  for  superconducting  magnetic  gradiometers  provided  surface 
speeds  exceed  1 rm/sec  or  so. 
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In  what  follows,  we  first  obtain  a general  solution  for 

electromagnetic  fields  induced  by  progressive  waves  in  a horizontally 

stratified  ocean.  We  then  restrict  consideration  to  wave  frequencies 

far  above  tidal  frequencies  and  present  solutions  for  transverse  electric 

type  fields  induced  by  surface  waves  and  by  internal  waves  in  an 

exponentially  stratified  ocean  and  an  ocean  stratified  in  two  layers. 

We  conclude  by  demonstrating  that  vertically  spaced  measurements 

of  noise  spectra  of  magnetic  field  gradients  above  an  ocean  surface 

offer  a unique  and  promising  means  of  obtaining  a measure  of  surface 

and  internal  wave  spectra  for  three  reasons:  (1)  proportionality  of 

magnetic  field  strengths  above  the  surface  to  a weighted  average  of 

sea  water  speed  over  an  ocean  depth  provides  a response  depending 

on  the  mode  structure  of  internal  waves,  (2)  exponential  decrease  of 

-kh 

field  strengths  with  height  above  the  surface  as  e provides  means 
of  wave  number  discrimination,  which  together  with  frequency  discri- 
mination offers  means  of  resolving  internal  wave  spectra  mode  by 
mode,  and  (3)  noise  spectra  of  magnetic  field  gradients  are  effectively 
cross-spectra  and  so  provide  directional  information  on  wave  spectra 
from  an  effective  point  measurement. 
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H.  ELECTROMAGNETIC  FIELD  GENERATED 


BY  MOVING  SEA  WATER 


Sea  water  moving  across  the  Earth’s  magnetic  field  induces 
an  electric  current  in  the  sea  that  generates  an  electromagnetic  field 
characteristic  of  the  velocity  field  in  the  sea.  Motion  of  sea  water 
with  velocity  V across  the  Earth’s  constant  magnetic  field,  B,  pro- 
duces a conduction  current,  7,  in  the  sea  that  we  express  in  MKS 
units  as 


J = <r(E  + VXB)  , (1) 

where  E is  the  electric  field  in  a stationary  frame  of  reference  and 
V X B is  the  electric  field  induced  by  sea  water  motion.  Electrical 
conductivity,  o-,  of  sea  water  is  about  4 mhos/m.  At  frequencies 
characteristic  of  sea  water  oscillations,  displacement  currents  in  the 
sea  are  negligibly  small  compared  to  conduction  currents  [Kraichman, 
1970],  so  that  we  determine  the  magnetic  field,  "b,  of  the  induced  cur- 
rent by  using  Ampere’s  law  in  the  form 

?Xb  = .uQJ  , (2) 

where  = in  X 10  henry/m.  A second  order  contribution  to  the 
induced  current  proportional  to  V X b is  negligible,  since  amplitudes 
of  induced  magnetic  fields  are  several  orders  of  magnitude  smaller 
than  the  Earth’s  field.  Finally,  electric  and  magnetic  fields  are 


6 


related  by  Faraday's  law,  so  that 

VXE  = - j|  . (3) 

Equations  (1'  (2),  and  (3)  together  with  appropriate  bound- 
ary conditions  determine  the  electromagnetic  field,  E and  b,  generated 
by  sea  water  motion  having  a velocity  field  V.  Because  electrical  con- 
ductivity of  sea  water  is  small,  the  Lorentz  body  force,  if  X B,  acting 
on  the  fluid  as  a result  on  the  induced  current  is  negligibly  small  com- 
pared to  pressure  and  buoyancy  forces.  Consequently,  hydromagnetic 
coupling  effects  are  absent. 

We  suppose  that  sea  water  flow  fields  are  described  by  motion 
of  an  incompressible  fluid  and  so  exclude  flow  fields  produced  by 
sound  waves.  The  divergence  of  velocity  then  vanishes,  and  the  vel- 
ocity field  is  the  curl  of  a vector  stream  function,  namely, 

V = VX*  , (4a) 

and  we  choose  the  gauge 

"V*  ^ = 0. 

The  electric  field  induced  by  moving  sea  water,  then,  is  the  sum, 

V X B = ?Xf$xB)  - v(^-B)  , (4b) 

of  transverse  and  longitudinal  type  fields. 
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A.  VECTOR  POTENTIAL  FORMULATION 

In  order  to  express  the  electromagnetic  field  in  terms  of  the 
vector  stream  function,  we  introduce  a vector  potential,  A,  for  the 
induced  magnetic  field,  so  that 

b = ?XA  . (5a) 

By  using  Faraday’s  law,  we  express  the  electric  field  as 

E = -|p  - vcf>  , (5b) 

where  <tis  a scalar  potential  function.  Because  the  divergence  of  J 
vanishes,  vector  and  scalar  potentials  are  related  by  the  expression 

V24>  + ■■■  - = -v2($.B)  . (5c) 

From  Ampere's  law,  we  obtain  the  additional  relation 

VXVXX  + = M°’VX( iXB) 

o o 

- nQar^(4>  + ^ • B)  (5d) 


between  potentials. 

By  specifying  the  gauge 

?•  A = -fiv  x , (6a) 


where 
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X 


(6b) 


= <t>  + B , 

we  obtain  the  expressions 

?2*  " VS?  = 0 (?a) 

and 

V2A  - Mo(r|^  = -V?X($XB)  (7b) 

from  Equations  (5c)  and  (5d),  respectively. 

Because  Equation  (7a)  is  homogeneous,  however,  we  make  a 
gauge  transformation  to  a gauge  in  which  V-  It  and  x vanish  and  express 
the  electric  field  and  current  in  terms  of  the  vector  stream  function 
and  a vector  potential  alone  as 

f = - |5  + v(*  • B)  (8a) 


and 


J = -o-~  + <r?X  (*X  B)  , (8b) 

where  the  vector  potential  A is  a solenoidal  solution  of  Equation  (7b) . 

Above  the  ocean  surface,  both  the  stream  function  and  elec- 
trical conductivity  vanish,  so  the  corresponding  vector  potential,  Ay, 
is  a solenoidal  solution  of  the  vector  form  of  Laplace's  equation. 
Below  the  ocean  floor,  the  stream  function  vanishes,  and  the  vector 

potential,  A , is  a solenoidal  solution  of  the  homogeneous  form  of 
s 
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Equation  (7b)  with  an  electrical  conductivity  or  corresponding  to  that 
of  suboceanic  strata.  The  electric  field  and  current  above  and  below 
the  ocean  are  given  by  corresponding  forms  of  Equations  (8a)  and  (8b), 
respectively.  Finally,  electric  and  magnetic  fields  produced  by  mov- 
ing sea  water  vanish  far  above  and  below  the  ocean,  and  their  horizontal 
components  are  continuous  across  the  surface  and  floor  of  the  ocean. 

B.  SCALAR  POTENTIAL  COMPONENTS 

In  order  to  determine  the  electromagnetic  field  generated  by 
a specified  velocity  field,  then,  we  determine  the  vector  potential  gen- 
erated by  the  corresponding  vector  stream  function,  so  that  derived 
electric  and  magnetic  fields  vanish  far  above  and  below  the  ocean  and 
have  continuous  horizontal  components.  Because  both  the  vector  stream 
function  and  vector  potential  are  solenoidal,  each  is  determined  by  two 
scalar  potential  components  (Morse  and  Feshbach.  1953].  For  our  pur- 
poses, we  represent  oceans  as  horizontally  stratified  and  so  express 
^and  A in  terms  of  scalar  potential  components  ^ ^ and  Aj,  Ag, 
respectively,  as 

= VX(^  + zX  V*2)  (9a) 

and 

A = V X (AjZ  + £XVA2)  , (9b) 


where  z is  a unit  vector  directed  vertically  downward. 

From  Equation  (7b),  we  then  find  that  each  pair  of  scalar 
potentials,  Aj,  and  Ag,  ^2,  is  related  by  the  expression 
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Furthermore,  each  scalar  component,  Ayj  and  A^g,  of  the  vector 
potential  above  the  ocean  surface  is  a solution  of  Laplace's  equation, 
and  each  scalar  component,  Agj  and  Agg,  of  the  vector  potential  below 
the  ocean  floor  is  a solution  of  the  homogeneous  form  of  Equation  (10). 
Because  scalar  potential  components  of  the  vector  stream  function 
and  vector  potential  are  paired,  each  component  of  the  vector  stream 
function  excites  its  own  mode  of  electromagnetic  field. 

C . ELECTROMAGNETIC  FIELD  GENERATED  BY  A 

PROGRESSIVE  WAVE 

In  order  tc  specify  the  mode  structure  of  electromagnetic 
fields  generated  by  ocean  waves,  we  express  each  of  the  scalar  poten- 
tial components  as  a Fourier  integral  of  a potential  function  corres- 
ponding to  a horizontally  progressing  ocean  wave  having  an  angular 
frequency  a>  and  horizontal  wave  vector!*.  For  example,  we  express 
scalar  potential  components  of  the  vector  stream  function  as- 

tfj(z,r,t)  = JJdwdI*^(z,i*,w)  e1^  , (11a) 

where  r is  a horizontal  coordinate  vector.  The  Fourier  coefficient 
^(z,k,o)  is  the  vertical  profile  of  a wave  progressing  horizontally  in 
a direction  1*  with  an  angular  frequency  w. 

Each  vector  field,  then,  is  a Fourier  integral  of  a corres- 
ponding vector  profile.  We  express  the  vector  stream  function,  for 
example,  as 


’W 


*(z,r,t)  =//d  dEf(z,Ti,«)  ei<"t‘5  r)  . (lib) 

By  using  Equations  (9a)  and  (11a),  we  express  the  vector  profile 
iA(z,  k,o>)in  terms  of  corresponding  potential  function  profiles  ^ and 

'l/2  as 

♦ = ik^n  + nXV^g)  . (11c) 


The  reduced  gradient  operator,  V , is  defined  by 

Vr  - z^z  - > (Hd) 

and  n = z X ic,  where  1c  is  a horizontal  unit  vector  pointing  in  the  direc- 
tion of  wave  propagation.  By  taking  the  curl  of  Equation  (lib),  then, 
we  obtain  the  expression 

~v  = ik(nV^2  - n X fj^'j)  (lie) 

for  the  vector  profile,  U(z,k,w),  of  the  velocity  field. 

Similarly,  vector  profiles,  a and  % of  the  vector  potential 
and  magnetic  field,  respectively,  are  expressed  as 

"a  = ik^n  + nXv^)  (12a) 


I an<* 

| = ik(nV^a2  - n X V^) 

j 


(12b) 
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Vector  profiles  of  electric  field  and  current,  T and  J,  are  given  by 
^ = cok  (a^n  + n XV^g)  + B)  (12c) 


and 


fij  = -ik[nV*aj  + n X %(^ag)]  . (12d) 

s 

Finally,  we  find  from  Equation  (10)  that  vertical  profiles  of 
scalar  components  of  the  vector  stream  function  and  vector  potential 
are  paired  according  to  the  relation 

2 — — 

V a.  - i/Mnua.  = -no-B’Vjt.  . (13) 

r j oj  0 r J 


Above  the  surface,  vertical  profiles,  ayl  and  av3,  are  solutions  of 
Laplace's  equation  with  a reduced  form  of  the  Laplacian  operator. 
Vertical  profiles  in  suboceanic  strata,  agj  and  agg,  are  solutions  of 
the  homogeneous  form  of  Equation  (13). 


As  is  evident  from  Equations  (12a),  (12b),  and  (12c),  the 
scalar  component  a j of  the  vector  potential,  and  hence  the  scalar 
component  ^ of  the  vector  stream  function,  generates  a transverse 
electric  type  field,  so  called  because  the  corresponding  electric  field 
vector  is  normal  to  the  plane  containing  the  magnetic  field  and  wave 
vectors.  The  scalar  component  ag  of  the  vector  potential,  and  hence 
the  scalar  component  ^ °*  the  vector  stream  function,  generates  a 
so-called  transverse  magnetic  type  field.  In  each  case,  the  electric 
field  also  contains  an  electrostatic  type  field  generated  by  the  gradient 
of  ^-5.  The  electrostatic  field  vector  lies  in  a vertical  plane  con- 
taining the  wave  vector. 
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We  observe  from  Equation  (12b)  that  the  transverse  magnetic 
mode  vanishes  above  the  surface,  because  scalar  components  of  the 
vector  potential  are  harmonic  above  the  surface.  Although  both  modes 
propagate  below  the  surface,  the  transverse  electric  mode  alone  propa- 
gates above  the  surface.  From  Equation  (12d),  we  also  note  that  the 
transverse  magnetic  mode  is  excited  by  electric  currents  in  a vertical 
plane  containing  the  wave  vector  and  that  the  transverse  electric  mode 
is  excited  by  horizontal  electric  currents  normal  to  the  plane. 

By  imposing  continuity  of  horizontal  components  of  electric 
and  magnetic  fields  at  the  ocean  surface  and  floor,  we  find  that 
modes  are  not  coupled  at  boundaries,  so  that  the  mode  structure 
of  electromagnetic  fields  generated  by  ocean  waves  is  determined 
by  the  vector  stream  function.  If  the  vector  stream  function  contains 
either  the  scalar  component  alone  or  the  scalar  component  ^ alone, 
then  the  electromagnetic  field  is  the  sum  of  a transverse  electric  type 
field  and  an  electrostatic  type  field,  in  the  first  case,  and  is  the  sum 
of  a transverse  magnetic  type  field  and  an  electrostatic  type  field,  in 
the  second  case.  In  each  case,  the  electrostatic  part  of  the  field 
results  from  surface  charge  required  to  halt  vertical  electric  currents 
at  the  surface.  For  vector  stream  functions  containing  both  scalar 
components,  the  electromagnetic  field  is  a sum  of  transverse  electric 
and  transverse  magnetic  type  fields  below  the  surface  and  a transverse 
electric  type  field  above  the  surface,  along  with  corresponding  electro- 
static type  fields. 

We  observe  from  Equation  (lie)  that  the  scalar  component 
'l' g of  the  vector  stream  function  produces  a velocity  component  nor- 
mal to  a vertical  plane  containing  the  wave  vector.  Consequently, 
velocity  fields  of  surface  and  internal  waves  at  frequencies  appreciably 
greater  than  tidal  frequencies,  for  which  the  velocity  vector  is  elliptically 
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polarized  in  a vertical  plane  containing  the  wave  vector  f Phillips,  1969], 
are  described  by  the  scalar  component  ^ alone.  Velocity  fields  of 
waves  with  frequencies  of  the  order  of  tidal  frequencies  contain  a 
transverse  component  of  velocity  and  so  generate  a transverse  mag- 
netic mode  as  well  as  a transverse  electric  mode.  The  transverse 
magnetic  mode,  however,  is  confined  below  the  surface. 

Because  we  are  interested  herein  in  internal  and  surface 
waves  with  frequencies  appreciably  greater  than  tidal  frequencies,  we 
confine  further  consideration  to  transverse  electric  type  fields. 


IH.  TRANSVERSE  ELECTRIC  TYPE  FIELDS  GENERATED 
BY  PROGRESSIVE  WAVES 


Progressive  waves  having  velocity  vectors  elliptically  polar- 
ized in  a vertical  plane  containing  the  wave  vector  generate  transverse 
electric  type  fields.  The  scalar  component  ^ of  016  vector  stream 
function  then  vanishes,  and  by  using  Equation  (lie),  we  express  the 
vector  profile  of  the  velocity  field  in  terms  of  circularly  polarized 
components  as 

v = ~ [(l+p)t)  + (l-p)0*]  , (14a) 

VI 

where  the  circularly  polarized  unit  vector  0 is  defined  by 


v = -p  (z  + i£)  , (14b) 

V2 

v*  is  the  complex  conjugate  of  v,  and  the  polarization  coefficient, 
p(z,  k,&»)  is  defined  by 


P 


1 d^z 
k dz 


(14c) 


The  vertical  component,  v _,  of  the  vector  profile  is  related  to  the 

z 

scalar  component  ^ of  the  vector  stream  function  by  the  expression 
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V 

z 


= k\  . (14d) 

The  unit  vector  u has  the  convenient  property  that  vX  n = iO.  Circu- 
lar polarization  of  v is  left-handed  with  respect  to  n,  and  that  of  v* 
is  right-handed. 

Vector  profiles  of  corresponding  transverse  electric  type 
fields  are  given  by  Equations  (12b)  and  (12c)  with  the  scalar  component 
a 2 of  the  vector  potential  profile  equal  to  zero.  The  remaining  scalar 
component,  a^,  of  the  vector  potential  profile  is  the  solution  of  Equa- 
tion (13)  that  gives  continuous  horizontal  components  of  electric  and 
magnetic  fields. 

Scalar  components  of  the  vector  potential  above  the  surface 
are  harmonic,  and  fields  vanish  far  above  the  surface.  Consequently, 
we  express  vector  profiles  of  electric  and  magnetic  fields  above  the 
surface,  t v and  0v,  in  terms  of  the  complex  amplitude  of  the  magnetic 
field  vector  profile  at  the  surface,  /3q,  as 

Kz 

«v  = *0¥> “P  n + ivZQ(B-n)  e*zV2v*  (15a) 

V La 

and 


/»„  - B0ekzC*  , (15b) 

where  z ^0  above  the  surface  and  i>ZQ  is  the  vertical  component  of  the 
velocity  profile  at  the  surface.  We  note  that  the  magnetic  field  above 
the  surface  is  circularly  polarized  with  a right-handed  sense  of  rota- 
tion with  respect  to  the  unit  vector  n. 


17 


Information  on  suboceanic  electrical  conductivity  profiles  is 
sparse.  Conductivity  of  the  Earth's  crust  below  the  oceans  is  much 
less  than  that  of  sea  water.  Conductivity  increases  with  depth  within 
the  Earth's  mantle  below  the  crust  and  becomes  much  greater  than 
that  of  sea  water.  Estimates  of  the  extent  of  the  poorly  conducting 
crustal  layer  beneath  the  oceans  range  from  a few  tens  of  kilometers 
to  a few  hundreds  of  kilometers  [Cox,  et  al. , 1970],  which  is  of  the 
order  of  ionospheric  heights.  For  simplicity,  then,  we  represent  the 
poorly  conducting  layer  as  a vacuum  and  consider  ocean  wavelengths 
that  are  much  less  than  the  extent  of  the  poorly  conducting  layer  and 
so  approximate  conductivity  profiles  of  suboceanic  strata  as  a vacuum 
region  of  indefinite  extent  below  the  ocean  floor. 


Consequently,  scalar  components  of  the  vector  potential  are 

harmonic  below  the  ocean  floor  as  well,  and  we  express  suboceanic 

vector  profiles  of  electric  and  magnetic  fields,  T and  /3  , in  terms  of 

s s 

the  complex  amplitude  of  the  magnetic  field  vector  profile  at  the 
bottom,  £d,  as 


e-k(z-D) 

vf 


A 

n 


(16a) 


and 


fl  e-k<z-D) 
;D  e 


A 

U 


(16b) 


where  z >D  and  D is  the  ocean  depth.  We  note  that  the  electrostatic 
part  of  the  electric  field  vanishes  in  suboceanic  strata,  because  the 
vertical  component  of  velocity  vanishes  at  the  ocean  floor. 

By  requiring  continuity  of  horizontal  components  of  electric 
and  magnetic  fields,  we  find  that  the  logarithmic  derivative  of  the 
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scalar  component  a^,  namely  (da^/dz)/a^,  equals  k at  the  ocean  sur- 
face and  -k  at  the  ocean  floor.  By  using  the  Green's  function 

g(M>  - -ie-k,2-{l  , »7a) 

we  write  Equation  (13)  for  the  scalar  component  a ^ as  the  integral 
equation 

aj({)  = «/B-Vr*j(ai)e'klz'cldz 
o 

*^T^/Dal^e"k,Z  ^kdz  ’ (17b) 

where 

* Mo- 

8=  “BT  * 

Expression  (17b)  is  a Fredholm  integral  equation  of  the  second  kind.  * 
We  write  the  solution  for  a^  as  the  Neuman  series  [Tricomi, 

1957] 


oo  m „ , , 

al  = 6Gf  + (“T>  8m+1Gm+1f  (18a) 

where 


kf  = B • V„  *1/1 
r l 


* The  trivial  solution,  aj=0,  is  the  only  solution  of  the  homogeneous 
form  of  Equation  (17b),  because  eigenvalues  of  a symmetric  kernel 
are  real  [Tricomi,  1957], 
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and  the  integral  operator  G is  defined  by 


(18b) 


The  quantity  5,  whose  reciprocal  is  a measure  of  the  speed  at  which 
fields  diffuse  over  a distance  of  one  ocean  wavelength,  is  given  in 
terms  of  ocean  wavelength,  X in  meters,  by 


8 


(4  X 10"7)A 


sec/m 


(18c) 


and  so  is  small  for  wavelengths  of  interest  herein.  Consequently,  we 
approximate  a^  to  first  order  in  8 by  retaining  only  the  first  term  in 
Equation  (18a),  which  is  an  integral  over  the  source. 

From  Equations  (12b),  (12c),  and  (12d),  we  then  obtain  the 
first  order  expressions 

?(£)  = -iS(|)n  J n(z)e_k^z“^  kdz  - (§-n)(Pxfi)  (19a) 

and 

/?({)  = -i 8v^[v/«.  n(z)e~k^~z^kdz 
o ’ 

+ v*  J e;  n(z)e”k^z"^kdz]  , (19b) 

for  vector  profiles  of  electric  and  magnetic  fields  in  the  ocean  (0££<D), 
and  the  first  order  expression 


3 


<rc. 
i)  n 


A 

n , 


(19c) 
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for  electric  current  in  the  ocean.  The  profile  of  the  transverse  compo- 
nent of  the  electric  field  induced  by  sea  water  motion,  « . (z),  is  given 

n 

by 


< = fi-PXB) 


(19d) 


From  Equations  (19a)  and  (19c),  we  note  that 


j - = r(fXB)  + iS(^)n  f €.  (z)e’’^z"^kdz 

k j i,  n 


(20) 


For  consistency,  then,  we  require  the  product  5(*>A)  to  be  small  as 
well.  Phase  speeds  of  ocean  waves  are  less  than  or  about  300  m/sec 
in  the  deepest  ocean,  where  D~  10  km,  so  that  a first  order  approxi- 
mation is  appropriate  for  ocean  wavelengths  less  than  several  kilo- 
meters, which  is  quite  adequate  for  our  purposes.  * 

We  conclude,  then,  that  the  electric  field  component  induced 
by  magnetic  field  oscillations,  which  is  normal  to  the  plane  of  polari- 
zation of  velocity,  does  not  appreciably  influence  electric  current,  so 
that  electric  current,  in  effect,  is  driven  by  the  transverse  component 
of  the  electric  field  induced  by  sea  water  motion.  Consequently, 
electric  current  and  magnetic  field  are  determined  directly  by 
the  velocity  field.  Magnitude  of  the  electric  current  is  proportional 
to  local  sea  water  speed.  Magnetic  field  strength  is  proportional  to 
a weighted  average  of  the  profile  of  sea  water  speed.  The  component 
of  the  electric  field  induced  by  sea  water  motion  that  is  in  the  plane 
of  polarization  of  velocity  simply  produces  surface  charges,  which 
halt  vertical  electric  currents. 


* Although  obtaining  exact  solutions  for  the  scalar  component  ai,  and 
hence  for  induced  fields,  is  straightforward,  resulting  expressions 
are  cumbersome  and  obscure. 
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A. 


VELOCITY  PROFILES 


In  a stratified  ocean  having  a Brunt-Vais%l&  frequency  pro- 
file, N(z),  the  vertical  component  of  the  vector  profile  of  velocity  of 
a progressive  wave  is  a solution  of  the  integral  equation 


N2(f)i/  «)G ({,  z)kd£ 

Zi 


v sinh  k(D-z) 
zo  sinh  kD 


(21a) 


which  is  equivalent  to  the  differential  formulation  given  by  Phillips, 
[1969].  The  symmetric  kernel  G(i , z)  is  given  by 


C%z) 


sinh  k£  sinh  k(D-z)  , _ 

SHhtB - ,or 


(21b) 


and  by 


G(£,z) 


sinh  kz  sinh  k(D-Q  ( 
sinh  kD  * 


(21c) 


The  horizontal  component  of  the  vector  profile  of  velocity,  i^(z),  is 
determined  by  the  relation 


1 

k dz 


(22a) 


since  the  divergence  of  velocity  vanishes.  Vertical  and  horizontal 
components  at  the  surface,  and  respectively,  are  related  by 
the  expression 


(£) 


l~l  V, 


ko 


(22b) 
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2 

where  g = 9.8  m/sec  , and  the  vertical  component  vanishes  at  the 
ocean  bottom.  The  dispersion  relation  for  progressive  waves  is  given 
by  the  expression 

2 < /*D9 

i^kc/tanh  kD  - cosh  kD  = J NlO^Osinh  k(D-£)kd£ . (23) 

Equation  (21a)  for  the  vertical  component  of  the  vector  pro- 
files of  velocity  of  a progressive  wave  is  a Fredholm  integral  equation 
of  the  second  kind,  so  that  properties  of  solutions  are  well  known 
[Tricomi,  1957],  Eigenfunctions  of  the  homogeneous  form  of  Equa- 
tion (21a)  correspond  to  internal  wave  modes  with  associated  eigen- 

o 

frequencies  u^.  At  internal  wave  frequencies,  the  ratio  w^/gk  is 
negligibly  small,  so  that  surface  displacements  effectively  vanish. 
Profiles  of  the  vertical  component  of  velocity  are  then  a superposition 
of  eigenfunctions  or  modes.  At  frequencies  much  greater  than  internal 
wave  frequencies,  the  profile  effectively  is  that  obtained  from  Equa- 
tion (21a)  for  a vanishing  Brunt- V&isal3.  frequency,  which  is  the  pro- 
file corresponding  to  surface  waves  on  a uniform  ocean. 

For  Brunt- VSisala  frequency  profiles  corresponding  to  actual 
ocean  stratification,  solutions  ordinarily  must  be  obtained  numerically 
or  by  approximate  methods  such  as  a WKB  approximation,  which  is 
appropriate  for  slowly  varying  profiles.  For  purposes  of  subsequent 
discussion  , however,  we  consider  three  model  profiles  that  repre- 
sent three  ideal  types  of  ocean  stratification:  (1)  a uniform  ocean, 

(2)  an  exponentially  stratified  ocean,  and  (3)  an  ocean  stratified  in 
two  layers. 

The  Brunt-Vaisala  frequency  vanishes  for  a uniform  ocean, 
and  surface  waves  are  then  the  only  mode  of  progressive  wave  motion. 
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The  corresponding  profile  of  the  vertical  component  of  velocity  and 
dispersion  relation  are  immediately  evident  from  Equations  (21a)  and 
(23),  respectively. 

The  frequency  profile  is  constant  or  uniform  for  an  exponen- 
tially stratified  ocean.  We  find  that  the  profile  of  the  vertical  compo- 
nent of  velocity  for  an  internal  wave  is  then  given  by  the  eigenfunction 
or  modal  expansion 

uz(z)  = \Q  n(-Dn  H?  sin  n*-(l-z/D)  (24a) 

and  that 


(24b) 


so  that  ^ko  n is  the  horizontal  component  of  velocity  at  the  surface 
resulting  from  the  nth  mode.  Corresponding  dispersion  relations  are 
given  by 


n = 1,  2,  3, . . . 


(24c) 


For  a two-layer  stratification,  we  represent  the  frequency 
profile  in  terms  of  a delta  function,  8(z-d),  centered  at  a depth  d,  so 
that 


N2(z)  = g(^W-d)  , (25a) 

where  LP/p  is  the  fractional  increase  in  density  between  layers,  which 

-3 

is  of  the  order  of  10  in  the  oceans.  By  neglecting  terms  of  the  order 
o 

of  (o  /gk,  we  then  obtain  from  Equation  (21a)  the  profile 
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uz(z)  = 1 % 


|iinh'k(D-^  ' G(d-Z) 
and,  from  Equation  (23),  the  dispersion  relation 


sinh  kD 


(25b) 


ai  = gk 


W[- 


sinh  kd  sinh  k(D- 
sinh  kD 


d)' 

“ ? 


(25c) 


for  a two-layer  stratification  representing  a sharp  thermocline  at  a 
depth  d.  We  note  that  the  velocity  profile  is  described  by  a single 
mode,  as  for  surface  waves  on  a uniform  ocean. 


B.  ELECTROMAGNETIC  FIELD  PROFILES 

The  profile  of  the  transverse  component  of  the  electric  field 
induced  by  sea  water  motion,  as  given  by  Equation  (19d),  is  propor- 
tional to  the  horizontal  component  of  velocity  in  polar  regions  and  to 
the  vertical  component  of  velocity  in  equatorial  regions.  As  a result, 
magnetic  fields  are  generated  by  the  horizontal  component  of  velocity 
in  polar  regions  and  by  the  vertical  component  of  velocity  in  equatorial 
regions.  At  intermediate  latitudes,  however,  both  components  of 
velocity  contribute,  and  fields  are  a superposition  of  fields  correspond- 
ing to  polar  and  equatorial  regions  with  relative  proportions  depending 
on  magnetic  latitude. 

Nonetheless,  the  horizontal  component  of  the  vector  profile 
of  velocity  is  proportional  to  the  vertical  gradient  of  the  vertical 
component,  because  the  divergence  of  velocity  vanishes,  so  that 
by  using  Equations  (14a)  and  (19b)  and  integration  by  parts,  we 
express  the  vector  profile  of  the  induced  magnetic  field  in  terms 
of  the  vertical  component,  u , alone  as 
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?(z)  = 8|v^(B.S)  |i/z(z)-vzoe"kj  -2(g.v)y,pz(0e“k(z"{)kd{|v 
-sj$B*z )uz(z)-2(B^)y*i.za)e"k(?“ZW  |>  , (26a) 


within  the  sea  (0  < z<  D),  and  as 


/3y(h)  = -Se 


-kh 


VSf(5-z)i/„„  - 2(B 


.D 


zo  _._>)/ ^({)e"k£kd{  t>  (26b) 


at  a height  h above  the  surface.  We  again  note  that  the  magnetic 

field  above  the  surface  is  circularly  polarized  and  that  its  magnitude 

~kh 

decreases  exponentially  with  height  as  e . Magnetic  field  strengths 

at  the  surfaco  depend  on  a weighted  integral  of  the  vertical  component 

of  velocity  over  an  ocean  depth  that  is  a Laplace  transform  of  the 

-k£ 

vertical  component  of  velocity.  Because  of  the  factor  e in  the 
integral,  the  contribution  of  deep  eddy  currents  to  magnetic  field 
strengths  at  the  surface  is  greater  at  smaller  wave  numbers.  Finally, 
we  observe  from  Equation  (26b)  that  the  gradient  in  a direction  u of  the 
vector  profile  of  the  magnetic  field  induced  above  the  surface  is  simply 
expressed  by 

(O-Vj.)  0V  » k >/2(fi>)  6V  . (26c) 

From  Equations  (19a)  and  (19b),  we  find  that  corresponding 
vector  profiles  of  electric  fields  are  given  by 

^ = “ *(5,n)l(i'*  *'*)»'  - (P-M *]  (26d) 

within  the  sea  and,  from  Equations  (15a)  and  (15b),  by 

% = + ^zo^-n)  "khvT0*  (26e) 
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above  the  surface.  We  note  that  the  impedance  of  the  transverse  elec- 
tric part  of  the  field  is  proportional  to  ocean  wave  phase  speeds.  The 
electrostatic  part  of  the  field  vanishes  in  polar  regions  and  is  circu- 
larly polarized  above  the  surface  with  a magnitude  proportional  to  the 
vertical  component  of  velocity  at  the  surface. 

Because  velocity  profiles  of  progressive  waves  in  a stratified 
ocean  depend  on  the  Brunt- V&isalS  frequency  profile,  description  of 
electromagnetic  fields  generated  by  progressive  waves  requires  a 
model  for  the  frequency  profile.  In  what  follows,  we  first  present 
results  for  electromagnetic  fields  generated  by  surface  waves  on  a 
uniform  ocean  and  then  for  electromagnetic  fields  generated  by  inter- 
nal waves  in  an  exponentially  stratified  ocean  and  an  ocean  stratified 
in  two  layers. 
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IV.  ELECTROMAGNETIC  FIELDS  GENERATED 
BY  SURFACE  WAVES 


By  evaluating  integrals  in  Equations  (26a)  and  (26b)  for  veloc- 
ity profiles  corresponding  to  surface  waves  on  a uniform  ocean  of 
depth  D,  we  find  that  vector  profiles  of  magnetic  fields  generated  by 
surface  waves  are  given  by 


P = 


-u  (TV 

o zo 
2k  sinh  kD 


j^B'y^e  k^sinh  kz  + (B*0)kzek^~z^]  0 


+ [(B-O)sinhk(D-z)  + (B-0*)k(D-z)e'k(D"z)]  0*j  (27a) 


within  the  sea  (0  < z < D)  and  by 


- -v 

ST  vzo 


[(B-0) 


sinh  kD 


(B-0*)] 


e’kh  0*  (27b) 


at  a height  h above  the  surface.  Corresponding  vector  profiles  of 
electric  fields  are  given  by 

7 ■ -s-Std  [-k(D'z)  +**-k(D-zi 


within  the  sea  and  by 


= ^(z,^y)n+  i^zo(B*n)e“kh  VT"0* 


(27d) 


above  the  surface. 
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We  observe  from  Equations  (27a)  and  (27c)  that  the  strength 
of  electric  and  magnetic  fields  generated  by  surface  waves  is  propor- 
tional to  the  vertical  component  of  sea  water  velocity  at  the  surface. 
For  small  wave  amplitudes,  then,  field  strengths  are  also  proportional 
to  wave  height,  £ t because  of  the  relation 

vZQ  = UQ  ^gk  tanh  kD  , (28a) 

which  expresses  the  kinematic  relation  between  vertical  displacement 
and  the  vertical  component  of  velocity  at  the  surface,  or  to  the  hori- 
zontal component  of  sea  water  velocity  at  the  surface,  because  of  the 
relation 


uzo  = ~l\o tanh  kD  ’ (28b) 

which  expresses  the  relation  between  velocity  components  at  the  surface . 

For  purposes  of  illustration,  we  choose  an  ocean  depth  of 

1 km  and  take  the  value  6. 24  X 10^  gamma  [Chapman,  1962]  for  the 

Earth’s  field  strength  in  polar  regions,  B , with  the  definition  one 
-9  2™ 

gamma  equals  10  webers/m  . We  express  the  Earth's  magnetic 
field  in  terms  of  the  unit  vectors  z,  k,  and  n as 

g - 

B = yjl  + 3 sin2 4 [z  sin  <t>  + (£  cos  n sin  6)  cos  <£]  , (29) 

where  $ is  the  angle  between  the  direction  of  wave  propagation  and 
magnetic  north,  4>  denotes  magnetic  latitude,  and  the  expression 
tan  <t>=  2 tan  * relates  dip  angle,  0 , and  magnetic  latitude.  Because 
wave  height  is  the  parameter  commonly  used  to  describe  surface  waves, 
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we  present  results  in  terms  of  field  strengths  generated  per  unit 
wave  height  and,  for  brevity,  consider  the  magnetic  field  alone. 

First,  however,  we  observe  from  Equation  (27b)  that  magni- 
tudes of  magnetic  fields  generated  at  the  surface  by  surface  waves  on 
a deep  ocean,  kD»l,  are  inversely  proportional  to  wave  number  for 
a given  surface  speed.  As  shown  by  the  solid  curves  in  Figure  1,  mag- 
netic field  strength  at  the  surface  increases  with  decreasing  wave  num- 
ber for  a constant  surface  speed,  because  the  contribution  of  deep  eddy 
currents  to  field  strength  at  the  surface  is  greater  at  smaller  wave 
numbers.  In  polar  regions,  electric  fields  induced  by  sea  water  motion 
are  proportional  to  the  horizontal  component  of  velocity,  and  magnetic 
field  strengths  at  the  surface  per  unit  horizontal  surface  speed  attain 
values  of  a few  gamma/(cm/sec)  at  wavelengths  much  greater  than  an 
ocean  depth.  In  equatorial  regions,  however,  electric  fields  induced 
by  sea  water  motion  are  proportional  to  the  vertical  component  of  vel- 
ocity, and  magnetic  field  strengths  at  the  surface  per  unit  vertical  sur- 
face speed  attain  values  of  several  tenths  of  a gamma/(cm/sec)  at  wave- 
lengths much  greater  than  an  ocean  depth.  Maximum  magnetic  field 
strengths  at  the  surface,  per  unit  surface  speed,  are  proportional  to 
ocean  depth  and  so  attain  values  of  the  order  of  10  gamma/(cm/sec) 
in  the  deepest  oceans,  for  which  D~  10  km. 

Horizontal  gradients  in  the  direction  of  wave  propagation  of 

the  magnetic  field  at  the  surface  (which  are  simply  expressed  as 

-ikjg,  however,  depend  only  on  the  product  of  wavenumber  and  ocean 

depth  for  a given  surface  speed.  As  shown  by  the  dashed  curves  in 

Figure  1,  magnitudes  of  horizontal  gradients  attain  maximum  values 
-3 

of  the  order  of  10  (gamma/m)/(cm/sec)  at  wavelengths  less  than  an 
ocean  depth  and  decrease  in  proportion  to  kD  at  wavelengths  much 
greater  than  an  ocean  depth. 
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Magnitudes  of  magnetic  fields  generated  at  the  surface  by  a 
one  meter  high  surface  wave  propagating  in  polar,  equatorial,  and 
intermediate  latitudes  (4>  = 90°,  0°,  and  45°)  are  shown  as  functions 
of  kD  by  the  solid  curves  in  Figure  2 and  correspond  to  the  left-hand 
coordinate  scale,  which  is  marked  in  units  of  gamma/m.  Dashed 
curves  delineate  corresponding  magnitudes  of  horizontal  gradients  in 
magnetic  field  strength  along  the  direction  of  wave  propagation  and 
correspond  to  the  right-hand  coordinate  scale,  which  is  marked  in 
units  of  (gamma/m)/m.  The  upper  abscissa  gives  the  ratio  of  ocean 
depth  to  wavelength,  A,  of  surface  waves.  Curves  corresponding  to 
equatorial  and  intermediate  latitudes  are  drawn  for  wave  propagation 
along  magnetic  meridians,  6=  0 or  n.  Induced  magnetic  fields  are 
independent  of  the  direction  of  wave  propagation  in  polar  regions  and 
are  proportional  to  cos  6 in  equatorial  regions.  At  intermediate  lati- 
tudes, the  dependence  on  direction  of  wave  propagation  is  only  signifi- 
cant in  a deep  ocean. 

As  is  evident  in  Figure  2,  magnetic  field  strengths  induced 
at  the  surface  by  a one  meter  high  surface  wave  increase  with  decreas- 
ing wave  number  and  attain  maximum  values  of  the  order  of  10  gamma/m 
in  polar  and  intermediate  latitudes  at  wavelengths  much  greater  than 
an  ocean  depth.  In  equatorial  regions,  however,  field  strengths 
decrease  with  decreasing  wave  number  for  wavelengths  much  greater 
than  an  ocean  depth,  because  the  vertical  component  of  sea  water 
velocity  at  the  surface  decreases  in  proportion  to  wave  number  for  a 
constant  wave  height,  as  is  evident  from  Equation  (28a).  Horizontal 
gradients  in  magnetic  field  strength  at  the  surface  increase  in  propor- 
tional to  VkD  for  wavelengths  less  than  an  ocean  depth  and  decrease  in 
proportion  to  kD  in  polar  regions  and  (kD)  in  equatorial  regions  for 
wavelengths  much  greater  than  an  ocean  depth.  Finally,  magnetic  field 
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strengths  per  unit  wave  height  shown  in  Figure  2 scale  according  to 
the  square  root  of  ocean  depth  in  kilometers,  but  magnitudes  of  hori- 
zontal gradients  per  unit  wave  height  seal''  according  to  the  inverse 
of  the  square  root  of  ocean  depth  in  kilometers. 

In  a deep  ocean  or  for  ocean  wavelengths  less  than  an  ocean 
depth,  we  find  from  Equations  (27a)  and  (28a)  that  vector  profiles  of 
magnetic  fields  generated  by  surface  waves  are  given  m terms  of  wave 
height  by  the  expression 

P = (5-k)  e~kz  (6*  + 2kz£)  (30a) 

within  the  sea  (z>0)  and  by  the  expression 

fiv  = ^ f(°)M  e"kh  (30b) 

at  a height  h above  the  surface.*  We  observe  that  variations  in  polari- 
zation and  magnetic  field  strength  with  depth  are  independent  of  latitude 
and  depend  only  on  the  ratio  of  depth  to  ocean  wavelength. 

Profiles  of  polarization  and  magnitude  of  the  magnetic  field 
induced  in  a deep  ocean  by  a one  meter  high  surface  wave  propagating 
in  polar  regions  with  a wavelength  of  10  m are  shown  in  Figure  3,  in 
which  altitude  and  depth  are  measured  in  units  of  an  ocean  wavelength. 
The  magnetic  field  vector  above  the  surface  is  circularly  polarized  in 
a right-handed  sense  with  respect  to  the  unit  vector  ft,  and  its  magni- 
tude decreases  exponentially  with  altitude.  Immediately  below  the  sur- 
face, the  magnetic  field  vector  is  elliptically  polarized  with  the  major 


* Expressions  for  a deep  ocean  approximation  agree  with  results 
presented  by  Weaver  [1965]. 
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axis  of  the  ellipse  aligned  vertically  and  with  a right-handed  sense  of 
rotation.  At  a depth  of  X/4ir,  the  field  is  linearly  polarized  in  a ver- 
tical direction.  At  greater  depths,  the  field  is  again  elliptically  polar- 
ized with  the  major  axis  of  the  ellipse  aligned  vertically,  but  the 
sense  of  rotation  is  left-handed.  At  great  depths,  polarization  becomes 
circular  with  a left-handed  sense  of  rotation. 

Maximum  magnetic  field  strengths  are  at  the  surface,  and 
values  scale  according  to  the  square  root  of  ocean  wavelength.  In 
equatorial  regions,  magnetic  field  strengths  are  less  than  field  strengths 
in  polar  regions  by  a factor  of  (1/2)  cos  0.  At  intermediate  latitudes, 
field  strengths  are  less  than  field  strengths  in  polar  regions  by  a fac- 
tor of 

l/2 

R = ^[(1  + 3 sin^<J>)  (sin^£  + cos^tf>  cos20)]  . 

In  shallow  seas  or  for  wavelengths  greater  than  an  ocean  depth, 
however,  variations  in  polarization  and  magnetic  field  strength  with 
depth  change  with  latitude.  Profiles  of  polarization  and  magnitude  of 
the  magnetic  field  induced  in  a shallow  sea  (kD  = 0. 1)  by  a one  meter 
high  surface  wave  propagating  in  polar  and  equatorial  regions  are  shown 
as  functions  of  z/D  in  Figure  4.  We  observe  that  profiles  are  symmetric 
in  polar  regions,  because  the  horizontal  component  of  sea  water  velocity 
is  effectively  independent  of  depth  in  a shallow  sea.  In  equatorial  regions, 
however,  profiles  are  asymmetric,  because  the  vertical  component  of 
sea  water  velocity  vanishes  at  the  ocean  floor. 
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V.  ELECTROMAGNETIC  FIELDS  GENERATED 
BY  INTERNAL  WAVES 


Description  of  electromagnetic  fields  generated  by  internal 
waves  is  more  complicated  than  for  surface  waves  because  internal 
waves  are  anistropic  and  propagate  at  oblique  angles  to  the  surface. 

As  a result,  numerous  internal  waves  or  modes  contribute  to  the 
velocity  field  corresponding  to  a particular  horizontal  wave  number. 
Electromagnetic  fields  generated  by  internal  waves,  then,  are  sums 
of  fields  generated  by  each  internal  wave  mode. 

Because  internal  waves  produce  negligibly  small  surface 
displacements,  we  neglect  the  vertical  component  of  velocity  at  the 
surface  and  find  from  Equations  (26a)  and  (26b)  that  vector  profiles 
of  magnetic  fields  generated  by  an  nth  mode  internal  wave  are  given 
by 

1Z  ) 

(B  • z)  uzn(z)-v^(B • v)£  VJl)emkizm°M  j v 

-V2(B.l)*)/°zn({)e'k({‘z)kd{Ji)»  , (32a) 


within  the  sea  (0  < z < D),  and  by 


Mo0-  -kh/s  a 


= “if"  e (B«P*)y*/  u__ (Oe  kd£  , 


vn 


-kt. 


(32b) 
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at  a height  h above  the  surface,  where  v(z)  is  the  vertical  component 
of  the  vector  profile  of  velocity  for  an  nth  mode  internal  wave.  Corre- 
sponding vector  profiles  of  electric  fields  generated  by  an  nth  mode 
internal  wave  are  given  by 

rn  - . <32c> 

within  the  sea,  and  by 


above  the  surface,  where  o»n(k)  is  the  eigenfrequency  of  an  nth  mode 
internal  wave  corresponding  to  a horizontal  wave  number  k.  We  note 
that  the  electrostatic  part  of  the  electric  field  vanishes  above  the  sur- 
face, since  the  vertical  component  of  velocity  ut  the  surface  is  neg- 
ligible. We  also  observe  that  the  magnetic  field  above  the  surface  is 
proportional  to  a Laplace  transform  of  the  vertical  component  of  the 
vector  profile  of  velocity. 

As  we  noted  previously,  velocity  profiles  corresponding  to 
actual  ocean  stratification  ordinarily  must  be  obtained  numerically 
or  by  suitable  approximate  methods.  Here,  we  consider  two  model 
profiles  that  represent  two  extremes  of  ocean  stratification:  (1)  an 
exponentially  stratified  ocean,  for  which  the  Brunt-VSis'dla  frequency 
profile  is  constant  and  (2)  an  ocean  stratified  in  two  layers,  for  which 
the  Brunt-VSisaia  frequency  profile  is  represented  by  a delta  function. 

A.  EXPONENTIAL  STRATIFICATION 

For  a constant  Brunt- Vaisala  frequency,  each  internal  wave 
mode  is  a sum  of  two  internal  waves  having  equal  amplitudes  and  prop- 
agating obliquely  to  the  surface  at  angles  #n  and  -#n,  where  tan  ^n=n^/kD, 
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with  a dispersion  given  by  the  relation  &>n/N=cos  0n,  which  is  equiva- 
lent to  Equation  (24c).  The  vector  profile  of  velocity  for  an  nth  mode 
internal  wave  is  then  expressed  in  terms  of  circularly  polarized  compo- 
nents as 


(-1)  V ( 

2iIL ) 0 sin  [mr(l-z/D)  + e ] 

V5T sin  6 ( 

n 


+ { Y sin  [n7r  (1-z/D)  - 0J  | , 


(33) 


where  is  the  horizontal  component  of  velocity  at  the  surface 
owing  to  the  nth  mode . 

Consequently,  we  find  from  Equations  (32a)  and  (32b)  that  the 
vector  profile  of  the  magnetic  field  induced  by  an  nth  mode  internal 
wave  is  given  by 


sin  (n7r(l-z/D)  + 2 0R)  - (-l)ne”kz  sin  26^ 


- (B*  0*)  sin  n^(l-z/D)|  0+  iAnj(B*0)  sin  nvil-z/D) 


(B-0*) 


sin  (mr(l-z/D)-20n)  + e"k(^"z) 


sin  20. 


n 


SO* 


within  the  sea  (0<  z<  D),  where 


(34b) 


(34a) 
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and  by 


'Vh)  ■ 


'Vn  IT 


1 - (-l)VkD 


1 + (nrr/kD) 


T 


e”kh  (B-v*)  v*  (34c) 


at  a height  h above  the  surface.  Corresponding  vector  profiles  of 
electric  fields  are  given  by  Equations  (32c)  and  (32d). 


For  purposes  of  illustration,  we  again  choose  an  ocean  depth 
of  1 km  and  consider  induced  magnetic  fields  alone.  We  present  results 
in  terms  of  magnetic  field  strengths  generated  per  unit  speed  at  the 
surface  on  a mode  by  mode  basis,  so  that  results  are  independent  of 
the  magnitude  of  the  Brunt- V3is515  frequency.  We  note,  however, 
that  surface  speed  and  the  maximum  vertical  displacement  of  an  nth 
mode  internal  wave,  £m  are  related  by  the  expression 


ko,n 


(-l)n+1N 


m,  n 


(35) 


4 


1 +(kD/mr)‘ 


which  results  from  the  kinematic  relation  between  vertical  displace- 
ment and  the  vertical  component  of  velocity.  Field  strengths  generated 
per  unit  maximum  displacement,  then,  are  proportional  to  the  Brunt - 
VSisam  frequency.  Furthermore,  high  order  modes  produce  greater 
surface  speeds  per  unit  max  imum  displacement  than  low  order  modes 

Magnitudes  of  magnetic  fields  and  corresponding  horizontal 
gradients  in  the  direction  of  wave  propagation  induced  at  the  surface 
by  internal  waves  having  surface  speeds  of  1 cm/sec  in  a 1 km  deep 
ocean  in  polar  regions  are  shown  as  functions  of  kD  by  solid  and  dashed 
curves,  respectively,  in  Figure  5 for  modes  1 and  2.  Magnitudes 
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of  horizontal  gradients  correspond  to  the  right-hand  coordinate  scale, 
marked  in  units  of  (gamma/m)/(cm/sec),  and  field  strengths,  to  the 
left-hand  coordinate  scale,  marked  in  units  of  gamma/(cm/sec). 

Field  strengths  generated  per  unit  surface  speed  shown  in  Figure  5 
scale  according  to  ocean  depth  in  kilometers,  but  magnitudes  of  hori- 
zontal gradients  depend  only  on  the  ratio  of  ocean  depth  to  wavelength, 
A.  At  equatorial  and  intermediate  latitudes,  field  strengths  at  the 
surface  are  less  than  field  strengths  in  polar  regions  by  the  factor  R 
as  given  by  Equation  (31). 

We  observe  from  Figure  5 and  Equation  (34c)  that  for  small 
values  of  n7r/kD,  magnetic  field  strengths  generated  per  unit  surface 
speed  are  inversely  proportional  to  wave  number  and  so  increase  with 
decreasing  wave  number,  because  deeper  eddy  currents  contribute  to 
field  strengths  at  the  surface  at  smaller  wave  numbers.  For  large 
values  of  nn/kD,  however,  magnetic  field  strengths  at  the  surface  are 
proportional  to  wave  number,  for  odd  numbered  modes,  and  to  the 
square  of  wave  number,  for  even  numbered  modes.  At  large  values 
of  n7r/kD,  sea  water  velocities  of  an  nth  mode  internal  wave  are 
nearly  linearly  polarized  and  horizontally  aligned  with  n phase  rever- 
sals over  an  ocean  depth.  In  polar  regions  then,  where  magnetic 
fields  are  generated  by  the  horizontal  component  of  sea  water  velocity, 
magnetic  field  strengths  generated  at  the  surface  decrease  with  wave 
number  because  contributions  of  deep  eddy  currents  increasingly 
oppose  contributions  of  shallow  eddy  currents.  In  equatorial  regions, 
however,  where  magnetic  fields  are  generated  by  the  vertical  compo- 
nent of  sea  water  velocity,  magnetic  field  strengths  generated  at  the 
surface  decrease  with  wave  number  both  because  contributions  of  deep 
eddy  currents  increasingly  oppose  contributions  of  shallow  eddy  cur- 
rents for  all  modes  except  the  fundamental  mode  and  because  the 
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vertical  component  of  velocity  decreases  with  wave  number  for  all 
modes  as  sea  water  velocity  approaches  linear  polarization. 

Profiles  of  magnitude  and  polarization  of  magnetic  fields 
induced  by  a fundamental  mode  Internal  wave  (n  = 1)  having  a surface 
speed  of  1 cm/sec  and  propagating  in  a 1 km  deep  ocean  in  polar 
regions  and  in  equatorial  regions,  along  magnetic  meridians,  are 
delineated  In  Figure  6,  for  kD  = 1.  We  observe  that  profiles  of 
field  magnitude  and  polarization  are  symmetrical,  as  are  velocity 
profiles.  The  major  axis  of  polarization  ellipses  for  odd  numbered 
modes  is  horizontal  in  polar  regions  and  vertical  in  equatorial 
regions.  For  even  numbered  modes,  however,  the  major  axis  is 
vertical  in  polar  regions  and  horizontal  in  equatorial  regions.  In 
each  case,  fields  are  linearly  polarized  at  a depth  of  D/2.  At 
intermediate  latitudes,  profiles  of  field  magnitude  and  polarization 
are  superpositions  of  profiles  corresponding  to  polar  and  equatorial 
regions,  with  relative  proportions  depending  on  latitude.  Nonethe- 
less, polarization  is  always  circular  both  above  the  ocean  surface 
and  below  the  ocean  floor,  and  field  magnitudes  decrease  expo- 
nentially as  e’k*1  with  height  above  the  surface  and  as  e“k(z_D) 
with  distance  below  the  floor. 

B.  TWO-LAYER  STRATIFICATION 

In  an  ocean  stratified  in  two  layers,  which  is  a model  com- 
monly used  to  represent  a sharp  thermocline,  only  a single  internal 
wave  mode  propagates.  Profiles  of  the  vertical  component  of  velocity 
corresponding  to  internal  waves  propagating  on  a sharp  thermocline 
are  expressed  by  Equation  (25b).  From  Equations  (32a)  and  (32b),  we 
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find  that  the  vector  profile  of  the  magnetic  field  induced  by  an  internal 
wave  propagating  on  a sharp  thermocline  at  a depth  d in  an  ocean  of 
depth  D is  expressed  by 
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above  the  thermocline  (0  < z < d),  where 


A ^ 

A1  * \o  "X"  ’ 


(36b) 


by 


/3g  = iAg  | (§•!)*)  sinh  k(D-z) 


+ (B-0) 


kze"kz  sinh  k(D-d)  - k(z-d)e“k^z"d^  sinh  kD 


sinh  kd 


A 

V 


iA, 


r 

(B-£*)k(D-z)e”k(D“z)  +(B-0)  sinhk(D-z) 


> , (36c) 


below  the  thermocline  (d  < z < D),  where 
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A _ * sinh  kd 
A2  " A1  sinh  k(Esd)  » 


(36d) 


and  by 


A 


= liv 


(w*\ 

(D/d  - l)e”k^  sinh  kd 

sinh  k(D-d) 

e~kh(B-i>*)i>  (36e) 


at  a height  h above  the  surface.  * Corresponding  vector  profiles  of 
electric  fields  are  expressed  by 


(37a) 


above  the  thermocline,  by 


e2  sinh  k(D-d) 

below  the  thermocline,  and  by 

?v  ■ I <H>®  , 


^ sinh  kd 
tA  ko 


+ e 


(B-n)r  k(D-z)  ^ 


(37b) 


(37c) 


at  a height  h above  the  surface,  with  a phase  speed  determined  by  the 
dispersion  relation  expressed  by  Equation  (25c). 

For  purposes  of  illustration,  we  choose  a thermocline  depth 
of  100  m and  an  ocean  depth  of  1 km  and,  as  before,  consider  induced 


* Expressions  for  the  vector  profile  of  the  magnetic  field  corresponding 
to  a deep  ocean  approximation  (kD»l)  agree  with  results  presented  by 
Beal  and  Weaver,  [1970], 
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magnetic  fields  alone.  We  present  results  in  terms  of  magnetic  field 
strengths  generated  per  unit  vertical  displacement,  of  the  thermo- 
cline.  Surface  speed  and  thermocline  displacement  are  related  by  the 
expression 


kd  sinh  k(D-d) 
sinh  kd  sinh  kf) 


(38) 


owing  to  the  kinematic  relation  between  vertical  displacement  and  the 
vertical  component  of  velocity.  We  choose  a representative  value  of 
0.001  for  the  fractional  increase  in  density  across  a thermocline. 

Magnitudes  of  magnetic  fields  and  horizontal  gradients  in  the 
direction  of  wave  propagation  generated  at  the  surface  by  internal  waves 
producing  aim  displacement  of  a sharp  thermocline  at  a depth  of 
100  m in  a 1 km  deep  ocean  in  polar  regions  are  shown  as  functions 
of  kD  in  Figure  7 by  solid  and  dashed  curves,  respectively.  For  a 
fixed  ratio  of  thermocline  depth  to  ocean  depth,  magnetic  field  strengths 
per  unit  thermocline  displacement  shown  in  Figure  7 scale  according 
to  the  square  root  of  thermocline  depth,  but  magnitudes  of  horizontal 
gradients  scale  according  to  the  inverse  square  root  of  thermocline 
depth.  At  equatorial  and  intermediate  latitudes,  magnitudes  of  mag- 
netic fields  and  field  gradients  at  the  surface  are  less  than  values  in 
polar  regions  by  the  factor  R expressed  by  Equation  (31). 

As  is  evident  from  Figure  7,  magnetic  field  strengths  gen- 
erated per  unit  thermocline  displacement  increase  with  decreasing 
wave  number  at  large  wave  numbers  (kD»l),  because  deeper  eddy 
currents  contribute  to  field  strengths  at  the  surface  at  smaller  wave 
numbers.  For  small  wave  numbers  (kD«l),  however,  field  strengths 
generated  at  the  surface  decrease  in  proportion  to  wave  number,  as 
is  evident  from  Equation  (36e) . At  small  wave  numbers,  sea  water 


42 


velocities  correspondint  to  an  internal  wave  propagating  on  a sharp 
thermocline  are  nearly  linearly  polarized  and  horizontally  aligned 
with  opposite  phases  above  and  below  the  thermocline.  In  polar  regions, 
then,  magnetic  field  strengths  generated  at  the  surface  decrease  with 
wave  number  because  contributions  of  eddy  currents  below  the  thermo- 
cline increasingly  oppose  contributions  of  eddy  currents  above  the 
thermocline  as  wave  number  decreases.  In  equatorial  regions,  how- 
ever, magnetic  field  strengths  generated  at  the  surface  decrease  with 
wave  number  because  the  vertical  component  of  velocity  decreases 
with  wave  number  as  sea  water  velocity  approaches  linear  polariza- 
tion. 

Profiles  of  magnetic  field  magnitude  and  polarization  in  polar 
and  equatorial  regions  are  shown  in  Figure  8,  for  wavelengths  much 
greater  than  a thermocline  depth  corresponding  to  kD  = 0. 5.  In 
polar  regions,  the  major  axis  of  polarization  ellipses  is  horizontal 
near  the  thermocline  and  becomes  vertical  away  from  the  thermo- 
cline. Three  polarization  reversals  occur  in  polar  regions:  one 
just  below  the  surface,  another  just  below  the  thermocline,  and 
the  last  just  above  the  bottom.  For  wavelengths  much  less  than  a 
thermocline  depth,  profiles  are  symmetrical  with  a polarization 
reversal  at  the  thermocline  and  one  directly  above  and,  below  the 
thermocline.  In  equatorial  regions,  the  major  axis  of  polarization 
ellipses  is  vertical,  and  a single  polarization  reversal  occurs 
directly  below  the  thermocline.  Profiles  are  again  symmetrical 
for  wavelengths  much  less  than  a thermocline  depth  with  the  polari- 
zation reversal  at  the  thermocline.  At  intermediate  latitudes,  profiles 
of  field  magnitude  and  polarization  again  are  superpositions  of  profiles 
corresponding  to  polar  and  equatorial  regions,  with  relative  proportions 
depending  on  latitude.  As  always,  polarization  is  circular  both  above 
and  below  an  ocean,  and  field  magnitudes  decrease  exponentially  as 
e ^ with  height  above  the  surface  and  as  e *^with  distance  below 
the  floor. 


43 


CONCLUSION 


As  we  have  demonstrated,  the  gradient  in  a direction  u 
of  the  magnetic  field,  la,  generated  at  a height  h above  an  ocean  sur- 
face by  a horizontally  progressing  ocean  is  expressed  by 

D 

(u-V)b  = -i/i0cre_kh  n*  (VxB)e"k^ kdC  , (39) 

— o 

where  V is  the  velocity  field  of  the  wave.  As  is  evident  from  Equa- 
tion (39),  gradients  of  magnetic  fields  generated  above  the  surface 
by  progressive  waves  are  circularly  polarized  in  a vertical  plane 
containing  the  direction  of  wave  propagation,  and  their  magnitudes 
decrease  exponentially  with  height  above  the  surface.  Magnitudes 
of  gradients  at  the  surface  are  proportional  to  a weighted  average  of 
sea  water  speed  over  an  ocean  depth. 

Because  magnitudes  of  magnetic  field  gradients  induced  at 
le  surface  are  proportional  to  a weighted  average  of  sea  water  speed 
rer  an  ocean  depth,  field  gradients  respond  to  the  mode  structure  of 
' "lal  waves,  as  is  also  evident  from  Figure  5.  Moreover,  because 
. aides  of  field  gradients  decrease  exponentially  with  height  above 
the  surface  as  e~^,  vertically  spaced  measurements  above  the  surface 
furnish  means  of  wave  number  discrimination,  which  together  with 
frequency  discrimination  affords  means  of  resolving  spectra  on  a 
mode  by  mode  basis.  Finally,  noise  spectra  of  magnetic  gradients, 
in  effect,  are  cross -spectra  and  so  provide  directional  information 
on  wave  spectra  from  an  effective  point  measurement.  The  imaginary 
part  of  the  spectrum,  for  example,  vanishes  for  isotropic  wave 
spectra.  Either  rotation  of  gradiometer  axes  or  multiple  sets  of  axes 
provide  information  on  anisotropy  of  a wave  field. 
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As  is  evident  from  Figures  1 and  5,  gradients  of  magnetic 

fields  generated  at  an  ocean  surface  by  both  surface  and  internal 

-3 

waves  are  of  the  order  of  10  gamma/m  for  surface  speeds  of 
1 cm/sec,  and  so  provide  sensible  signals  for  recently  developed 
superconducting,  magnetic  gradiometers  at  surface  speeds  typically 
produced  by  ocean  waves.  Magnetic  field  gradients  produced  at  an 
ocean  surface  by  fluctuating  ionospheric  currents  are  of  the  order  of 
10"4  gamma/m  or  less,  so  that  magnitudes  of  gradients  generated 
by  ocean  waves  exceed  noise  levels  produced  by  ionospheric  currents 
by  a factor  of  ten  or  more  for  typical  surface  speeds.  Measurements 
of  fluctuating  magnetic  field  gradients  at  a fixed  point  above  an  ocean 
surface,  then,  provide  an  estimate  of  ocean  wave  spectra. 

Nonetheless,  a hallmark  of  oceanographic  observations 
is  the  absence  of  a stationary  platform.  Use  of  a so-called  transverse 
axis  gradiometer  (in  which  the  two  pickup  loops  lie  in  one  plane)  to 
record  fluctuating  gradients  of  field  components  transverse  to  Earth's 
magnetic  field,  however,  suppresses  apparent  gradient  fluctuations 
owing  to  changes  in  gradiometer  orientation.  By  aligning  the  gradio- 
meter axis  so  that  a line  of  force  of  Earth's  field  passes  through  the 
center  of  each  pickup  loop,  changes  in  gradiometer  orientation  then 

produce  apparent  gradient  fluctuations  that  are  of  the  order  of  the 

-2 

product  of  Earth's  magnetic  field  gradient  (-10  gamma/m)  and  the 
sine  of  the  angular  deviation  in  alignment,  or  about  10"^  gamma/m 
for  a deviation  of  one  milliradian. 

Thus,  we  submit  that  using  appropriately  designed,  supercon- 
ducting, magnetic  gradiometers  to  make  vertically  spaced  measure- 
ments of  magnetic  field  gradients  generated  above  the  surface  by  ocean 
waves  provides  a unique  and  promising  means  of  obtaining  a measure 
of  ocean  wave  spectra  --  for  both  surface  and  internal  waves. 
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FIGURE  CAPTIONS 


Figure  1.  Magnitudes  of  magnetic  fields  and  horizontal  gradients  in 
the  direction  of  wave  propagation  generated  at  the  surface  of  a 1 km 
deep  ocean  by  surface  waves  having  horizontal  surface  speeds  of 
1 cm/sec  in  polar  regions  and  vertical  surface  speeds  of  1 cm/sec 
in  equatorial  regions  and  propagating  along  magnetic  meridians  in 
equatorial  regions  shown  as  functions  of  kD  along  the  lower  abscissa 
and  as  functions  of  D/X  along  the  upper  abscissa.  Solid  curves 
delineate  magnetic  field  strengths  per  unit  surface  speed  and  corres- 
pond to  the  left-hand  coordinate  scale  marked  in  units  of  gamma/(cm/ 
sec).  Dashed  curves  delineate  magnitudes  of  horizontal  gradients 
along  the  direction  of  wave  propagation  per  unit  surface  speed  and 
correspond  to  the  right-hand  coordinate  scale  marked  in  units  of 
(gamma/m)/(cm/sec). 

Figure  2.  Magnitudes  of  magnetic  fields  and  horizontal  gradients  in 
the  direction  of  wave  propagation  generated  at  the  surface  of  a 1 km 
deep  ocean  by  a one  meter  high  surface  wave  propagating  in  polar 
regions  ( $ = 90°)  and  along  magnetic  meridians  at  intermediate  and 
equatorial  latitudes  (<t>=  45°  and  0°)  shown  as  functions  of  kD  along 
the  lower  abscissa  and  as  functions  of  D/X  along  the  upper  abscissa. 
Solid  curves  delineate  magnetic  field  strengths  per  unit  wave  height 
and  correspond  to  the  left-hand  coordinate  scale  marked  in  units  of 
gamma/m.  Dashed  curves  delineate  magnitudes  of  horizontal 
gradients  along  the  direction  of  wave  propagation  per  unit  wave 
height  and  correspond  to  the  right-hand  coordinate  scale  marked  in 
units  of  (gamma/m)/m. 
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Figure  3.  Profiles  of  polarization  and  magnitude  of  the  magnetic  field 
generated  in  a deep  ocean  (kD»l)  by  a one  meter  high  surface  wave 
having  a wavelength  of  10m  and  propagating  in  polar  regions.  Altitude 
and  depth  are  measured  in  units  of  an  ocean  wavelength. 

Figure  4.  Profiles  of  polarization  and  magnitude  of  the  magnetic 
field  generated  in  a 1 km  deep  ocean  by  a one  meter  high  surface 
wave  propagating  in  polar  regions  and  along  magnetic  meridians  in 
equatorial  regions  and  having  a wavelength  corresponding  to  kD  = 0. 1. 
Altitude  and  depth  are  measured  in  units  of  an  ocean  depth,  D. 

Figure  5.  Magnitudes  of  magnetic  fields  and  corresponding  horizontal 
gradients  in  the  direction  of  wave  propagation  generated  at  the 
surface  of  a 1 km  deep  ocean  in  polar  regions  by  internal  wave 
modes  1 and  2 each  having  surface  speeds  of  1 cm/sec,  shown  as 
functions  of  kD  along  the  lower  abscissa  and  as  functions  of  D/x 
along  the  upper  abscissa.  Solid  curves  delineate  magnetic  field 
strengths  per  unit  surface  speed  and  correspond  to  the  left-hand 
coordinate  scale  marked  in  units  of  gamma/(cm/sec).  Dashed  curves 
delineate  magnitudes  of  horizontal  gradients  in  the  direction  of  wave 
propagation  and  correspond  to  the  right-hand  coordinate  scale  marked 
in  units  of  (gamma/m )/(c m/sec). 

Figure  6.  Profiles  of  polarization  and  magnitude  of  the  magnetic 
field  generated  in  a 1 km  deep  ocean  by  a fundamental  mode  internal 

wave  (n=l)  propagating  in  polar  regions  and  along  magnetic  meridians 
in  equatorial  regions  and  having  a surface  speed  of  1 cm/sec  and  a 
wavelength  corresponding  to  kD  = 1.  Altitude  and  depth  are  measured 
in  units  of  an  ocean  depth. 


51 


Figure  7.  Magnitudes  of  'magnetic  fields  and  horizontal  gradients  in 
the  direction  of  wave  propagation  generated  at  the  surface  of  a 1 km 
deep  ocean  by  internal  waves  producing  aim  displacement  of  a sharp 
thermocline  at  a depth  of  100  m in  polar  regions  shown  as  functions  of 
kD  along  the  lower  abscissa  and  as  functions  of  D/x  along  the  upper 
abscissa.  Solid  curves  delineate  magnetic  field  strengths  per  unit 
displacement  and  correspond  to  the  left-hand  coordinate  scale  marked 
in  units  of  gamma/m.  Dashed  curves  delineate  magnitudes  of  hori- 
zontal gradients  in  the  direction  of  wave  propagation  and  correspond 
to  the  right-hand  coordinate  scale  marked  in  units  of  (gamma/m)/m. 

Figure  8.  Profiles  of  polarization  and  magnitude  of  the  magnetic  field 
generated  in  a 1 km  deep  ocean  by  an  internal  wave  propagating  in  polar 
regions  and  along  magnetic  meridians  in  equatorial  regions  and  pro- 
ducing aim  displacement  of  a sharp  thermocline  at  a depth  of  100  m 
at  a wavelength  corresponding  to  kD  = 0. 5.  Altitude  and  depth  are 
measured  in  units  of  a thermocline  depth. 
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